Frontiers in Neuroimaging

Multi-echo functional magnetic resonance imaging (fMRI) acquires signals with distinct contrast profiles across echo times (TEs). Empirical evidence suggests that different TEs capture distinct signal contributions with varying sensitivity to BOLD and non-BOLD processes from long to short TEs (Stroman, 2002; Krishnamurthy, 2023; Dong, 2024). Previous literature has not yet examined how brain states extracted from temporal dynamic functional connectivity (dFNC) evolve across TEs. Here we present a methodological and empirical study of TE-dependent evolution of dynamic connectivity states, with a phenomenological two-component model for interpretation. Multiple echoes were acquired with an echo planar time-resolved imaging (EPTI) sequence. TE-specific components were extracted via atlas-based group ICA (Neuromark 1.0), and brain connectivity states were generated via windowed dFNC and spatial ICA clustering of windows. Within-subject and group-level states evolved across TEs in network organization and variance. States exhibited TE-dependent variance profiles, with distinct peak echo times; and group-level points of state transition were independently mapped to TE-dependent transitions in variance. We modified a bi-exponential model to simulate and model changes in the dFNC time-course across TEs. The bi-exponential model revealed high consistency with subject states (R2 fit=0.86, res. error=0.55, 86% model convergence rate), suggesting a two-component signal system contributing to dynamic TE-dependent connectivity changes in state. A significant correlation between relative model terms was estimated on a cell-wise basis across states (r(24) = -0.43, p=0.03, 95% CI [-0.70, -0.04]). This suggests that specific bi-exponential signal parameters (i.e., short FNC1,TE=0, long FNC2,TE=0, and long FNC2,decay) were associated with the observed TE-dependent variations in connectivity states across TEs. In sum, our analyses provide downstream implications for TE-aware state analysis dependent on connectivity magnitude (FNCTE=0) and rate of connectivity decay (FNCdecay), in a two-component signal model, with future implications for TE-dependent state biomarkers.

Background: Ultra-low field (ULF) magnetic resonance imaging (MRI) offers accessible point-of-care neuroimaging, yet its capability to support multivariate network analyses remains understudied. Methods: We systematically evaluated the feasibility and reliability of Morphometric Similarity Networks (MSN) and Morphometric Inverse Divergence (MIND) networks derived from paired high-field (HF, 3T/1.5T) and ULF (64mT) scans in a primary cohort(N=23) and an independent replication cohort (N=19). We tested the impact of super-resolution processing, cortical parcellation granularity (Desikan-Killiany vs. Destrieux), and mesh resolution (fsnative vs. fsaverage5) on edge-level cross-scanner agreement (Lin's concordance correlation coefficient, CCC) and subject-level identifiability. Results: MIND networks demonstrated higher edge-level HF-ULF agreement (mean CCC = 0.95) and ULF reproducibility (mean CCC = 0.98) compared to MSNs (mean CCC = 0.86-0.87 and 0.91-0.94, respectively). In contrast, MSNs yielded higher subject-level differential identifiability across field strengths due to a broader distinction between within- and between-subject similarity. Reconstructing MIND networks on native high-resolution meshes (fsnative) maximised cross-scanner agreement, whereas coarser meshes and finer cortical parcellations enhanced individual network identifiability. Conclusions: Structural similarity networks can be reliably constructed from ULF MRI data. Optimal analytic choices depend on study objectives: MIND on native meshes maximises cross-scanner stability, whereas MSNs excel at preserving individual uniqueness, paving the way for globally scalable connectomic research.

Quantitative magnetic resonance imaging (qMRI) supports both basic neuroscience and clinical research, providing the foundation for longitudinal studies of brain development, aging, and neurological disease. However, its validity rests on a little-studied assumption: that physical relocation of imaging hardware leaves underlying neuroimaging biomarkers intact. When research facilities expand or modernize, an entire scanner undergoes radical disruption--field ramp-down, physical transport, shielding reconstruction, and re-shimming. Whether qMRI metrics remain stable following such invasive procedures without introducing systematic biases is largely unquantified. Here, we systematically evaluated the impact of relocating a 3T clinical scanner over a distance of 30 km on a broad, multimodal suite of neuroimaging markers in healthy volunteers. To disentangle variance components, we benchmarked within- and cross-site stability using percent variation (PV; measurement precision), intraclass correlation coefficients (ICC; absolute reliability), and Spearman rank correlations (R2; spatial and order agreement). We show that scanner relocation introduces no greater variance than baseline within-site assessments. High-resolution brain structural morphometry from different sequences (multi-echo MPRAGE, compressed sensing MPR2RAGE) and quantitative relaxometry (compressed sensing MP2RAGE T1 mapping) demonstrated excellent stability (PV<1%, ICC=0.89-0.99, R2>0.93). Diffusion tensor (mean diffusivity: MD, fractional anisotropy: FA) and mean signal kurtosis (MSK) metrics from multi-shell diffusion MRI maintained low variability across white and deep gray matter (global cross-site ICC>0.90; global PV: MD<1.7%, MSK<1.2%, FA<3.8%; global R2: MSK>0.93; MD>0.90; FA>0.77), confirming relocation robustness for diffusion-derived microstructural markers. Resting-state functional networks displayed high cross-site reproducibility (ICC>0.90), while task-based fMRI activation peaks retained acceptable spatial and magnitude agreement (ICC=0.30-0.60) despite intrinsic physiological variability. In the absence of software or hardware upgrades, physical scanner relocation leaves qMRI biomarkers remarkably intact--providing the empirical foundation required to preserve data integrity across major facility transitions and multi-site longitudinal cohorts. Key points that summarize the paperO_LIHardware transitions do not compromise qMRI data integrity: Physical relocation of a 3T scanner (including field ramp-down, 30 km transport, and re-shimming) introduces no greater variance in quantitative MRI metrics than baseline within-site test-retest assessments. C_LIO_LIMultimodal biomarker robustness: Brain morphometry, T1 relaxometry, diffusion tensor and kurtosis, resting-state functional connectivity display comparable within- and cross-site stability reflected by different metrics (percent variation, intraclass correlation coefficient and Spearman correlation) C_LIO_LIEmpirical backing for longitudinal research: In the absence of hardware or software upgrades, multimodal qMRI biomarkers remain remarkably stable, reassuring researchers that long-term cohort studies and multi-site harmonization can survive major facility relocations. C_LI Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=175 SRC="FIGDIR/small/26363342v1_ufig1.gif" ALT="Figure 1"> View larger version (57K): org.highwire.dtl.DTLVardef@773d03org.highwire.dtl.DTLVardef@46a09borg.highwire.dtl.DTLVardef@da5bc0org.highwire.dtl.DTLVardef@1bdbf9d_HPS_FORMAT_FIGEXP M_FIG This study evaluates the impact of physically relocating a 3T MRI scanner on neuroimaging biomarkers. By comparing before and after relocation sessions on healthy volunteers, we conclude that the procedure does not introduce significant variability. Parameters for morphometry, T1 relaxometry, diffusion, and resting-state functional connectivity show high stability. Task-based fMRI showed greater within-site variability, but no increase in variability following relocation, consistent with the inherently higher variability of task-evoked responses.. The results confirm the integrity of qMRI data after infrastructure relocation, ensuring reliability for longitudinal studies. C_FIG

The choroid plexus (ChP) is a highly vascularized structure located in the brain ventricles that plays important roles in multiple central nervous system functions. ChP enlargement has been implicated in many neuropsychiatric conditions such as multiple sclerosis, depression, schizophrenia, and Alzheimer's disease. However, limitation in the accuracy of the current segmentation tools highlight the need for more reliable methods to quantify ChP volume in healthy and clinical human populations. The shape and size of the ChP vary drastically from person to person, and currently there are no large ChP datasets which allow comprehensive training of segmentation models. Here we introduce DataPlex, a collection of 700 structural MRI brain images where the lateral ventricle ChP has been manually labeled and checked by expert human raters. DataPlex is balanced and representative across a wide range of scanning parameters, demographics, and clinical profiles, and subject-level information will be available for each labeled image. DataPlex will be open source to support the creation of reliable segmentation tools for the human lateral ventricle ChP.

There is increasing demand for quantitative measurements in neuroimaging research, including absolute quantification of regional myelin density. Fast whole-brain macromolecular proton fraction (MPF) imaging enables acquisition and reconstruction of quantitative maps of MPF, a biomarker of myelin, and has been used in many research applications. To date, no standard protocols have been published for parcellating MPF maps. FreeSurfer is a widely used software package for automated brain parcellation. It is designed for use with high-resolution T1-weighted anatomical images, typically acquired with an MPRAGE sequence. MPF maps inherently provide high contrast between white and gray matter and can potentially be used as input images for parcellation, enabling simultaneous quantitative tissue characterization and volumetric assessment. However, it is unclear how the use of MPF maps as input images with FreeSurfer software affect reliability and repeatability of regional MPF and volume estimates. We evaluated three FreeSurfer parcellation workflows for this purpose. The workflows differed in the amount of neuroimaging data required and computational intensity. MPF and MPRAGE data were acquired in two separate sessions for 11 adults. Reliability was evaluated using mean relative differences in estimates, Bland-Altman analysis, and intra-class correlation coefficients (ICCs). Repeatability was assessed using within-subject coefficients of variation (CVws). All three workflows produced similarly high within-subject repeatability across scans. However, reliability in regional MPF estimates was lower in gray matter than in white matter for workflows that used MPF maps as input. Workflows that use MPF maps as FreeSurfer input may be adequate for estimating mean MPF in GM and WM of cortical parcels and in subcortical regions. However, the MPRAGE-based workflow is recommended when reliable individual-level estimates or volume estimates are required.

Background and purpose: Accurate identification of cardiac MRI volume orientations is essential for reliable image interpretation and for enabling downstream automated analysis pipelines. However, orientation labels and commonly assigned manually, making the process time-consuming and prone to variability. Recent advances in vision-language models offers new opportunities for automated label generation by jointly modeling visual and semantic information. Methods: In this study, we systematically compare a conventional end-to-end transformer-based captioning model with a CLIPCap architecture for automated cardiac MRI view labelling. Both models were evaluated on a dataset of 1184 de-identified cardiac MR images spanning across six clinically relevant orientations: short axis, two-, three-, four-chamber views, left and right ventricular outflow tract views. The caption transformer was trained from scratch, whereas CLIPCap leveraged pretrained CLIP image embeddings combined with lightweight mapping network and a frozen language model. Results: Experimental results demonstrate a substantial performance advantage from CLIPCap, which achieved an overall accuracy of 98% compared to 68% for the caption transformer. CLIPCap consistently delivered high precision and recall across all orientation classes, including those of high clinical importance, while the end to end transformer exhibited unstable performance and complete failure in certain views. Conclusion: These findings highlight the benefits of leveraging large scale multimodal pretrained model for medical image labeling tasks, particularly in data limited settings. The study suggests that CLIPCap provides a more robust and clinically reliable solution for automated cardiac MRI orientation labeling, supporting its integration into efficient and scalable clinical workflows.

PurposeWith growing interest in pediatric MRS to understand the neurochemistry of brain development, it is essential to establish pediatric-specific metabolite T1 and T2 relaxation times for reliable quantification of metabolites. Furthermore, characterization of age-related changes in relaxation times may provide insight into developmental differences in tissue microstructure and cellular composition. Therefore, the goals of this study were two-fold: (1) to estimate metabolite T1 and T2 relaxation times in a pediatric cohort; and (2) to assess linear relationships between age and relaxation times. MethodsTI- and TE-series data were acquired from the bilateral thalamus in 20 and 25 participants aged between 0 and 7 years, respectively. Data were processed and modelled in Osprey 2.5.0. Metabolite T1 and T2 relaxation times of total N-acetylaspartate at 2.0 ppm (tNAA2.0), total creatine at 3.0 ppm (tCr3.0), total creatine at 3.9 ppm (tCr3.9), total choline (tCho), myo-inositol (mI) and the combined glutamate and glutamine (Glx) were estimated by fitting the modelled signal amplitudes to TI-series inversion recovery signal and TE-series exponential decay equations, respectively. Linear relationships with age were tested for metabolite T1 and T2 relaxation times using Spearman correlations. ResultstNAA2.0 and Glx showed significant changes in T1 with age. tNAA2.0, tCr3.0, tCho and mI showed significant decrease in T2 with age. ConclusionSeveral metabolites showed significant changes in relaxation times with age in the thalamus region. Therefore, use of age-specific relaxation correction for pediatric studies using the linear equations provided is encouraged for reliable quantification of metabolites.

Neuroimaging research focusing on subcortical brain regions showed that they play a role in a wide range of cognitive functions and associated their structural alterations to various psychiatric disorders. Methods to analyse such structures were developed to observe changes, not simply in terms of volumes, but in terms of surfaces, which are able to detect subtle and local changes in shape and geometrical configuration. Such advanced analyses have however been largely restricted to command line environments with relatively limited access to Python users or researchers with low technical expertise. We present SubCortexMesh as a user-friendly toolbox which covers automated surface estimation from popular subcortical volume segmentations (FreeSurfer and the Functional Magnetic Resonance Imaging of the Brain Software Library (FSL)), computes shape-related vertex-wise metrics (thickness, surface area and curvature) for whole cohorts and includes statistical analyses. This is all done inside Python, automatically computing all subjects of a given preprocessed directory, with the explicit intent to minimise steps and manual coding required from the users. SubCortexMesh includes statistical tools which allow conventional random field theory-based cluster analyses on native subject metrics, standardised within common surface templates, thus containing a whole workflow necessary to run an up-to-date shape-wise subcortical analyses in one package.

Introduction: This study aimed to characterize the neurochemical profile of Alzheimer s disease (AD) continuum stages using proton magnetic resonance spectroscopy (1H-MRS) data for brain metabolites and macromolecules. We used Orthogonal Projection to Latent Structures Discriminant Analysis (OPLS-DA) to classify diagnostic groups and assess feature importance. Methods: Our dataset consisted of 7 T short-echo time 1H-MRS spectra of the posterior cingulate cortex/precuneus of individuals diagnosed with AD (n = 27) or Mild Cognitive Impairment (MCI, n = 27) and Healthy Controls (HC, n = 31). Classification used 100-round Monte Carlo cross-validation based on metabolite concentrations and macromolecular peak intensities. Variable importance in projection (VIP) was assessed across classification models to identify the most relevant neurochemical features. The resulting metabolite rankings were compared with previous literature to evaluate their biological plausibility. Results: Classification was strongest for AD versus HC (AUC=0.88) and lower for comparisons involving MCI. Metabolites contributing most strongly to classification included myo-inositol, ascorbate, and choline-containing compounds, consistent with previous studies on AD-related neurochemical alterations. In addition, specific macromolecular components were identified as relevant features, suggesting a potential contribution of MM signals to the characterization of AD-related neurochemistry. Discussion: These findings support the use of short-echo-time 1H-MRS as a complementary tool for characterizing neurochemical changes in the AD continuum, with potential applications in disease characterization and monitoring, warranting validation in larger cohorts.

Background: Whole-brain segmentation of neonatal low-field MRI is needed for regional volumetric analysis, but established segmentation systems are designed around different field strengths, contrasts, and age ranges. Methods: We first applied complete end-to-end pipelines, including FreeSurfer, FSL-based processing, Infant FreeSurfer, FastSurfer, BrainSuite, SPM12/CAT12, iBEAT v2, the dHCP structural pipeline, and M-CRIB-S, directly to 0.35T neonatal T2-weighted images. Because those systems did not produce a complete 47-region output, we next tested released end-to-end AI segmentation models, including SynthSeg, FastSurfer, FastSurferVINN, and QuickNAT. These were highly sensitive to input format, contrast, voxel geometry, and the high-field training domain. We then decomposed the task into skull stripping, registration, super-resolution, atlas propagation, and tissue refinement and tested the available components individually. Results: No complete end-to-end pipeline or released AI model delivered a usable 47-region segmentation in direct application. Stage-wise testing showed that preprocessing, reconstruction, registration, atlas propagation, and expectation-maximization each solved only part of the problem. Their final assembly produced complete 47-region outputs across the neonatal cohort. Conclusions: Direct application of established pipelines was insufficient for 0.35T neonatal T2-weighted MRI. A stage-wise combination of compatible components provided the workable solution.

The hippocampus is particularly vulnerable to hypoxia, yet freedivers, who repeatedly undergo voluntary breath-hold hypoxia, do not appear to show the memory impairments commonly reported after involuntary hypoxic exposure. We investigated whether freediving training is associated with functional connectivity reorganization at the whole-brain and hippocampal levels, and whether hippocampal connectivity relates to preserved episodic memory performance. Seventeen male recreational freedivers underwent fMRI before and after a 7-month standardized training period, and 20 matched aerobic-trained controls were assessed at baseline. Functional connectivity analyses examined (i) whole-brain ROI-to-ROI connectivity, (ii) seed-to-voxel connectivity from the bilateral hippocampus and selected cortical seeds, and (iii) associations between hippocampal connectivity and episodic memory. Compared with controls, freedivers showed distinct connectivity patterns during apnea and rest. After training, whole-brain analyses revealed reorganization involving salience, frontoparietal, sensorimotor, visual, and cerebellar systems. Hippocampal analyses showed reduced connectivity with sensorimotor and default mode regions, alongside increased connectivity with cerebellar and visual areas, with stronger effects after training. In freedivers, specific post-training hippocampal connectivity was associated with episodic memory performance. These findings suggest that repeated voluntary hypoxia during freediving training is associated with selective functional reorganization of hippocampal and large-scale brain networks. This pattern may reflect adaptive neuroplasticity linked to preserved episodic memory under intermittent hypoxic exposure.

Automatic segmentation of medical images could accelerate clinical and research workflows, yet existing tools are rarely compared independently on shared benchmarks, making it difficult to assess genuine progress. This is particularly consequential for limb muscle segmentation in magnetic resonance imaging (MRI), which underpins volumetric analysis, radiomics, and fat fraction quantification used as biomarkers in neuromuscular disease -- applications where segmentation errors can directly affect interpretation. Despite the availability of multiple open-source segmentation tools spanning diverse architectures, no independent head-to-head comparison has appeared in the literature. We address this gap by evaluating eight tools, six of which run independently, and the other two of which are meant to be auxiliary algorithms once muscles have been somehow initially labeled. We evaluated these tools on three-dimensional MRI volumes across three base datasets and then later against derived datasets, assessing both quantitative metrics and qualitative usability. Our results reveal substantial performance variation across methods and datasets, with many methods showing reduced accuracy on pathological cases. Domain-specific models trained on large datasets consistently outperformed foundation models and newer general-purpose architectures. Our results reveal practical trade-offs between accuracy, generalizability, and ease of use that are critical for clinical adoption, while raising concerns for deployment on underrepresented populations and pathologies.

Background: Interoception refers to the processing, integration, interpretation, and regulation of bodily signals by the brain. It is crucial to motivational and affective functioning such as subjective experience of arousal (i.e., intensity of emotion). Previous studies found altered self- report of arousal in people with schizophrenia, but the underlying biological mechanism is unclear. This study aimed to investigate interoceptive processing and emotional experience in schizophrenia by triangulating self-report, cardiac dynamics, and heartbeat-evoked potential (HEP; an event-related potential that reflects cortical processing of heartbeats). Methods: A total of 39 participants with schizophrenia spectrum disorders (SZ; 36% female) and 31 healthy controls (HC; 35% female) reported subjective arousal level when viewing images of varying emotional intensity while electroencephalogram and electrocardiogram were being recorded. Results: We found that intense images led to larger heartrate decelerations and higher subjective arousal ratings in HC, but did not modulate heartrate changes in SZ and were less predictive of their subjective ratings. Subjective ratings were predictive of cardiac dynamics in HC, but not in SZ. SZ showed a less negative HEP amplitude than HC when viewing less intense images. Lastly, SZ with less severe positive and negative symptoms exhibited cardiac dynamics more similar to those of HC. Conclusions: These findings suggest both reduced cardiac reactivity to emotional stimuli and a disconnect between interoceptive afferent signals and cortical processing, leading to different subjective emotional experience in people with schizophrenia. Studying interoception using a multimodal approach may reveal important insights into the underlying mechanisms of altered emotional experience in schizophrenia.

Resting-state connectivity can predict task-evoked fMRI activation, but correspondence with an individual task map may partly reflect a shared population pattern. We evaluated the Variational Resting-state-to-Task Prediction TransformeR (VRPTR), a three-dimensional encoder-decoder combining a compressed Transformer bottleneck, variational latent sampling, and multiscale skip connections. Models were trained in 360 healthy adults from the WU-Minn Human Connectome Project and evaluated in 40 held-out participants from the same dataset for the story-versus-math language contrast. VRPTR achieved mean voxel-map Pearson r=0.642 and Dice AUC=0.519, exceeding compact volumetric BrainSurfCNN-like and SWIFUN-like comparators by {Delta}r=0.0376 and 0.0335, respectively. A template fixed from the 360 training maps yielded residual correlation 0.364; VRPTR improved correlation but not mean absolute or root mean squared error over this template. A matched factorial ablation associated Transformer processing with improved map correlation ({Delta}r=0.0168); the variational main effect did not survive multiple-comparison correction. Component swaps indicated that subject-specific activation patterns depended more on multiscale skip features than on the bottleneck representation. Latent sampling produced a modest association between posterior standard deviation and voxelwise absolute error (mean Spearman {rho}=0.091). Raw 95% intervals covered only 4.7% of observed values. Five-fold calibration within the held-out cohort increased coverage to 94.8%, and spatially varying uncertainty modestly improved probabilistic scores over a constant-width model. Uncertainty also increased under input degradation. These findings support the contribution of Transformer processing to language-map prediction and distinguish it from subject-specific spatial conditioning and relative uncertainty estimation. Independent validation of calibrated uncertainty and evaluation in clinical cohorts remain necessary.

Identifying progressive dementia during its prodromal stage is a challenging yet critical task for healthcare providers. Traditional clinical assessments are often not sensitive enough to detect warning signs of cognitive decline. Machine learning (ML) methods have shown promise in modeling clinical trajectories associated with neurodegenerative disorders such as Alzheimer's disease (AD). Validation of these models' forecasting abilities, however, is widely lacking in the literature. To address this gap, we established a longitudinal cohort from the National Alzheimer's Coordinating Center (n = 15,727) comprising individuals who were cognitively normal (CN) or had mild cognitive impairment (MCI) at baseline. Cognitive status and clinical observations collected over up to 10 visits were used to train several ML classifiers. We employed the gradient-boosted tree ensemble method XGBoost with a custom feature-engineering pipeline. On non-imputed data, our models achieved an AUC-ROC of 0.9447 [0.9417-0.9476] and PR-AUC of 0.8554 [0.8497-0.8611] for the CN to MCI task, and an AUC-ROC of 0.9059 [0.8983-0.9134] and PR-AUC of 0.9214 [0.9156-0.9272] for the MCI to AD task. In a lead-time forecasting analysis spanning mean follow-up periods of 56.1 months (range, 9.1-214.1) for CN to MCI and 42.7 months (range, 9.1-126.1) for MCI to AD, the non-imputed models detected future progression before diagnosis in 29.58% [28.29%-30.87%] of CN to MCI progressors and 22.41% [21.36%-23.46%] of MCI to AD progressors. Among these individuals, the mean lead time was 33.04 months [32.52-33.56] and 17.95 months [17.71-18.19], respectively. False alarm rates among non-progressors were 12.6% [11.5%-13.7%] and 20.0% [18.4%-21.6%], respectively. These findings suggest that longitudinal machine learning models can identify future cognitive decline before diagnosis, indicating clinical utility for early dementia risk stratification.

We conducted an open-label pilot study of bilateral nucleus accumbens-anterior limb of the internal capsule (NAc-ALIC) deep brain stimulation (DBS) in three patients with treatment-resistant schizophrenia (TRS). Clinical symptoms, global functioning, cognition, stimulation parameters and intracranial local field potential (LFP) were assessed over 24 weeks. Two patients met the prespecified clinical response criterion, defined as a [&ge;]40% reduction in Positive and Negative Syndrome Scale (PANSS) total score, by week 8 and showed PANSS reductions of 69.0% and 69.4%, respectively, at week 24. The third patient showed a 18.7% reduction at week 24. Responders showed improvements in global functioning and selected cognitive domains. Longitudinal LFP recordings revealed increased gamma-band relative power and reduced theta-band activity in responders, whereas these changes were less consistent in the non-responder. DBS was generally well tolerated, with transient stimulation-related adverse effects resolving after reprogramming. These findings provide preliminary clinical evidence supporting further investigation of NAc-ALIC DBS in TRS and suggest that individualized programming and intracranial physiological measures may contribute to the identification of treatment-response signatures.

Submillimeter functional magnetic resonance imaging (fMRI) based on blood-oxygenation-level-dependent (BOLD) signal enables the study of brain function at the submillimeter level, uncovering insights into fine-scale organisations like cortical layers and columns. However, its inherently low contrast-to-noise ratio (CNR) and signal-to-noise ratio (SNR) often limit its reliability and applicability. Noise Reduction with Distribution Corrected Principal Components Analysis (NORDIC PCA) is a locally low-rank denoising algorithm that reduces thermal noise levels in BOLD fMRI in a local patch manner. However, local patches often contain a mixture of signals from multiple tissues that negatively affects the low-rank structure of the patches, which limits the denoising capabilities of the algorithm. We propose an alternative approach for patch formation by gathering similar non-local voxels, dubbed voxel-matching (VM) NORDIC. The results on submillimeter-resolution BOLD fMRI data indicate that VM-NORDIC effectively promotes the low rank of the patches by boosting signal redundancy, allowing for more efficient noise attenuation. Moreover, the method barely affects spatial smoothness due to the non-local voxel selection based on time-series similarity. In particular, VM-NORDIC outperforms standard NORDIC with default local patching (Standard-NORDIC) in terms of temporal SNR (tSNR) (~9-90% larger than Standard-NORDIC; ~23-250% larger than the original) and spatial smoothness estimates (~20% of the smoothness induced by Standard-NORDIC). These improvements are fundamental to improving the validity and precision of fMRI studies at submillimeter resolutions.

Background and Hypothesis: First-line antipsychotics are ineffective in ~40% of people with first episode psychosis (FEP). Electroencephalogram (EEG)-based predictive biomarkers could be a scalable solution to reduce trial-and-error. Study Design: We acquired EEG in 56 minimally treated FEP individuals before treatment initiation, and 34 matched healthy controls (HC). Treatment response was defined as >50% improvement in BPRS scores by 12 weeks of treatment. Modified k-means clustering of resting-state EEG yielded six microstates including canonical topographies (A-D). The first principal component (PC1) of microstate properties derived from HCs served as a normative reference onto which FEP patients were projected. Baseline PC1 was used to predict treatment response in adjusted regression analyses. To test stability of results, an independent dataset of HCs (n=92) was fitted using the same microstate solution and compared to our native HC dataset. PC1 scores derived from this independent normative dataset were also used for predicting response in FEP. Study Results: Individual microstates replicated established patient-control differences (e.g., higher presence of microstate C compared to D, p<0.001). PC1 explained 34.8% of variance. Baseline PC1 predicted subsequent treatment response (adjusted OR=0.63, p<0.001), yielding AUC=72.3% in leave-one-out cross-validation (pperm<0.001) when combined with clinical information. PC1 from our HC transported well to the second healthy control EEG dataset (n=92) (r=0.9,p<0.001; r=0.89,p<0.001 ). PC1 in FEP derived using PC1 weights from a second normative dataset also predicted response in FEP (OR=0.70, p=0.03). Conclusions: EEG microstate architecture differing from the expected norm is associated with lower odds of treatment response in FEP.

Researchers increasingly have access to neuroimaging data with repeated measures within subjects. Linear mixed-effects modeling offers a valuable way to analyze such data, but its application remains limited in the neuroimaging literature due to high computational cost and a lack of mainstream analysis packages. Here we benchmark a recently published algorithm, the Fast and Efficient Mixed-Effects Algorithm (FEMA), against a conventional implementation (MATLAB's fitlme, restricted maximum likelihood) across 94,830 functional connectivity edges. In contrast to the Adolescent Brain Cognitive Development (ABCD) Study Dataset used in the development and testing of FEMA, the dataset used in the current analysis has a smaller sample size (n = 378) and comprises individuals with chronic pelvic pain scanned at up to four visits over three years. Despite these differences, the two implementations produced highly concordant results: fixed-effect estimates and test statistics correlated at r [&ge;] 0.999 (Lin's concordance correlation coefficient [&ge;] 0.999), and the two approaches reached the same statistical conclusion for 99.98% of edges. FEMA completed each analysis roughly 20 times faster. Agreement was weaker for the estimated variance components, where FEMA attributed systematically less variance to the participant random effect, yielding a slightly lower intraclass correlation in approximately 82% of edges. We further show that most residual disagreement reflected FEMA's default variance-grid resolution rather than the estimator itself, and could be eliminated at negligible computational cost. These findings support the use of FEMA for connectome-wide analysis of repeated-measures neuroimaging data.

Background: Traditional neuropathological assessment of Alzheimer's disease (AD) relies on manual, semi-quantitative scoring of amyloid-beta (A{beta}) burden. This approach is limited by high inter-rater variability and often qualitative, failing to capture subtle, continuous pathology changes. While 11C-PiB PET neuroimaging provides non-invasive longitudinal monitoring, it has a limited dynamic range and lacks the sub-millimeter spatial resolution necessary to track localized microvascular and parenchymal A{beta} trajectories. Methods: We introduce ViTamyloid, a generalizable modular computational pathology pipeline for high-throughput A{beta} segmentation and classification in postmortem tissue. In Stage 1, a pathologist-supervised hybrid framework combining an EfficientNet-B4 encoder with a Residual U-Net decoder, complement with a localization-aware Faster R-CNN detector, to generate binary segmentation masks and isolated lesion candidates. In Stage 2, classification model powered by a frozen, high-resolution UNI2-h computational pathology foundation backbone and a ViT-Adapter, classifies candidates into a fine-grained four-class taxonomy: classic cored plaques (CCP), coarse-grained plaques (CGP), capillary cerebral amyloid angiopathy (CAA Type 1), and parenchymal/leptomeningeal CAA (CAA Type 2). We bench-marked out-of-distribution (OOD) robustness on external datasets from different brain banks, and validated clinical scalability across 1,646 whole-slide images from three brain banks. Results: The ViTamyloid classifier achieved strong feature-separation metrics, yielding area under the receiver operating characteristic curve (AUROC) and precision-recall curve (AUPRC) values of 0.994/0.958 for CCP, 0.989/0.945 for CGP, 0.992/0.977 for CAA Type 1, and 0.997/0.978 for CAA Type 2. More importantly, ViTamyloid consistently outperforms other pipelines' models and frozen features in all external datasets with different antibodies, and scanner settings, and image quality. Conclusion: By replacing empirical, from-scratch architectures with a generalizable foundation model adapter framework, ViTamyloid provides a highly sensitive scalable digital pathomics toolkit. This pipeline expose continuous microscopic neuropathological dynamics beyond the spatial and sensitivity plateaus of in vivo neuroimaging, offering objective, quantitative endpoints for Alzheimer's disease clinical trials and disease subtyping.